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O ' I. INTRODUCTION 



Some physical phenomenon happen to be expressible from topological properties of specific models. The integer 
^ Hall conductivity is one of them |^ , [ p7[ . In a simple model of non- interacting electrons moving in a two-dimensional 
periodic potential V{x,y) subject to a uniform perpendicular magnetic field B^ and a low electric field Ey, the Hall 
conductivity Cxy of a given filled Landau electronic band turns out to be proportional to an integer C: 

, C is the Chern index of the band, describing the topology of its fiber bundle structure [||. 
' For a better understanding of this phenomenon, and to bring out the conditions of possible experimental measures, 
^ I we investigate in this paper the value of C as a function of the potential V. This is done by semi-classical methods, 
I ' and the tunnelling effect appears to be the fundamental mechanism for a non-zero conductivity. 
' O [ In the limit of high magnetic field B^ , the above model is mapped onto the well-known Harper model, by means of 
^ the Peierls substitution: the potential V is considered as a perturbation of the cyclotron motion, and the averaging 
method of mechanics gives an effective Hamiltonian equal to the average of V on the cyclotron circles. We neglect 
the coupling between the Landau bands |^^. For a high magnetic field, (hence for a small cyclotron radius,) this 
transformation gives an effective Hamiltonian Heff{q,p) ~ V{q,p), which is bipcriodic in position and momentum 
(the phase space is a 2D torus), and an effective Planck constant /le// = hc/{eB). In this approximation, trajectories 
5J] ' are the levels lines of Heff ~ V. Furthermore, the expression of he/f shows that the high magnetic field regime 
Ci ' corresponds to the semi-classical limit. This model will be the starting point of our study in the next section. 

Quantum mechanics on the torus has been extensively used as a convenient framework to study basic properties of the 
semi-classical limit like quantum chaos ( @, [|l^, [D), or the tunnelling effect ( |31 3^, [D). For this purpose, a 
convenient tool is the Bargmann and Husimi representation which maps a quantum state to a function on the phase 
space. These representations are constructed with coherent states and will be recalled in section 3. 
The new results presented in this paper are the conditions under which the tunnelling effect between three different 
wells can be responsible for a non-zero Hall conductivity. The conditions will be expressed by specifying the special 
positions the three wells must have inside a periodic cell. 

It is worth mentioning that due to its topological aspect, it is natural to study the Chern index values in a generic 
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situation, because topological properties are robust against perturbations. Secondly, in the generic ensemble of 
Hamiltonian under study, the integer values of the Chern indices are characterized by the position of the boundary 
where their value changes by one unit. These boundaries turn out to correspond to degeneracies in the spectrum 
We are therefore brought to study the generic location of these degeneracies, in the space of Hamiltonians. This is 
done in section 5 and 6. 

Under reasonable assumptions on the tunnelling interaction in phase space, we find that the boundary of constant 
Chern index domain form ellipses in any generic two-dimensional subspace of the Hamiltonian's space. These analytical 
results are well confirmed by numerical calculations in section 7. 

These results extend previous work by the first author for the Hall conductivity resulting from the tunneling effect 
between two wells in a given periodic cell ( (lO|, (l^). Although the methods looks similar, calculations and results 
are quite different. 

II. QUANTUM MECHANICS ON THE TORUS 

A. Classical dynamics on the torus 

Let us consider a one-degree of freedom Hamiltonian (hence an integrable dynamics) , periodic both in position q 
and momentum p, with respective periods Q and P : 

H{q,p) = H{q + Q,p) = H{q,p + P) (1) 

The function H{q,p) can be decomposed into its Fourier series: 

H{q,p) ^ ^ c„i,„2 exp(i27rni^)exp(i27rn2^) (2) 

This decomposition will be used for numerical computations, but is not essential for the theoritical part of this work. 
Since _ff is a real valued function, the complex coefficients c„j„2 must satisfy: 

Crii,ri2 = c-,ii -„2 eff*, {ni,n2) e 'Z'^ 

The trajectories q{t),p{t) evolve on the plane, but because of periodicity, they can be considered as trajectories on 
the Torus Tqp, obtained by identifying the opposite sides of the cell [0, Q] x [0, P]. 
For example, the trajectories of the Harper Hamiltonian 

H{q,p) = - cos(2^^) - i cos(2^|) (3) 

are displayed in figure ol p. 0. Some energy levels are contractibles, others are not. 




Q ; 2Q ^ 



FIG. 1. Trajectories of the Harper model (^. 
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B. Quantum mechanics on the torus 



The Hilbert space of a particle with one-degree of freedom is (M) , with the fundamental operators of position q 
and momentum p acting on a wave function s-lso noted \ip > ■ 

We now choose a symmetric quantization procedure: to the classical Hamiltonian H, we associate the quantum 
Hamiltonian: 

H= -c„j^ „2 exp(i27rni^) exp(i27rn2-^) + hermitian conjugate 

We denote by Tq [respectively Tp] the translation operator by one period. Tq translates by Q a wave function 
ipix) and Tp translates its Fourier transform 'ijj{p) by P: 

fgijix) = ij{x - Q) (4) 
fpi^ip) = i^ip - P) (5) 

We may rewrite equations (Q) and ^ as: 

Tq = exp{-iQp/h), Tp = exp{iPq/h) 
Quantum mechanically speaking, the periodicity Eq. reads: 

[H,fQ]^[H,fp]^0 (6) 

To continue, we now have to assume that 

[fQ,fp]=0. (7) 

It is easy to prove that: 

TqTp = e-'^fpTQ 
hence (0) is fulfilled if and only if there exists an integer N such that: 

A-^geW (8) 

This hypothesis can be regarded as a geometric quantization condition, which states that there is an integer 
number of Planck cells in the phase space. The semi-classical limit h ^ corresponds to the limit N +oo. 

Remark 1 For a given periodic Hamiltonian H and a given Planck constant h, the ratio is generically not an 
integer. But the spectral properties of H are in some sense continuous with respect to h, hence we may approximate h 
by h so that |^ becomes a rational ^ ( /fJ^j. Now to fulfill our hypothesis, we consider H as a periodic hamiltonian 
with periods {Q, Dx P) (or equivalently [D x Q, P), many other combinations are possible if D is not a prime number). 
In the sequel, P and Q denotes periods of H and not necessarily primitive periods of H. From this, we see that it is 
natural to expect additional translation symmetries of the Hamiltonian inside the ( non primitive ) periodicity cell and 
the main results of this work will apply to the D = 3 case. 

We now assume that Hypothesis (||) is fulfilled. 

According to the commutation relations (H) and (0) , the Hilbert space {M) may be decomposed as a direct sum 
of the eigenspaces of the translation operators Tq and Tp: 

L''{]R)= J J HNiei,d2)d0id92 (9) 

HNiOuO.) = ||V. > such that I >- > 
I I Tp\ip exp(i6'2)|V' > 
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with {91,62) G [0, 27rp related to the periodicity of the wave function under translations by an elementary cell. The 
space of the parameters {01,62) has also the topology of a torus, and will be denoted by Tg. 

The space Hn{6i,62) is not a subspace of L^{IR), the space of physical states, it is a space of distributions in 
the X representation, we will sec later that it is a subspace of a weighted space in the Bargmann representation. 
We will now show that i/jv(^'i, ^'2) is finite dimensional. Let \i} >G Hn{6i,62). The Fourier transform of ip{x) is 
02-Floquet-periodic of period P, so ip{^) is discrete, it is a sum of Dirac distributions supported at points distant 
from = from each other. Moreover, '>p{x) is 6'1-periodic, hence -0 is characterized by the N coefficients at the N 
Dirac distributions supporting points in the interval q G [0, Q[. Eventually we get: 

dimo- Hn{6i,62) - N 

Because of Eq. the Hamiltonian H is block-diagonal with respect to the decomposition Eq. (^), so we have to 
consider the spectrum a{6i, 62) of H. The operator H acts on Hn{6i, 62) a.s a N x N hermitian matrix, its spectrum 
is made of eigenvalues: 

a{6i,62) = {Ei{6i,62),...,En{6i,62)} 

Let ■(/'jv be the corresponding eigenfunctions: 

H\M&uS2) >= En{6i,62)\M&i,02) > n=l,..,N 

For a given level n, as {61,62) G Tg are varying, the energy level £^{61,62) form a band, and (assuming that £^{61,62) 
is never degenerate \/6), the eigenvectors \^n{6i,62) > form a 2D surface in the quantum states space. But for a fixed 
{61,62), and any A e (T, AjV'nl^i, 6*2) > is also an eigenvector. So the family of eigenvectors for the level n form a 
complex-line-bundle (of fibre = (T 9 A ) in the projective space of the bundle — > Tg. The topology of this line 
bundle is characterized by an integer C„ S Z, called the Chern index. Because of the natural Hilbert scalar product 
on L'^{1R) = J J Hm{6i,62) ddid62, which induce the Berry (or Chern) connection ( (2^), this topological number 
is explicitly given by the integral of the Berry (or Chern) curvature ( p7|]): 

^" " / deM0i,02)\de,i^^{6u62) >-< deM0i,02)\deM6i,62) >) d6id62 (10) 

This expression has been used intensively for our numerical calculations. Moreover, it can be shown (see e.g. Q, |^) 
that: 

N 

E^" = i (11) 

n=l 



III. HUSIMI AND BARGMANN REPRESENTATION 



We have seen previously that the space Hn is not a subspace of L^{1R). For functions belonging to Hm, it will be 
more useful to introduce a phase space representation of the quantum states, called the Bargmann representation ( 

Consider a Quantum state IV' >£ L^{1R). In order to characterize the localization of \ip > in the phase space near 
the point {q,p), we first construct a Gaussian wave packet \qp > (coherent state) defined in the x-representation by: 



/ 1 i {x~ q)^ 
< x\qp ( j exp(-pa;) exp( — — ) 

The notation \qp > recalls that the coherent state is localized (in the semi-classical limit) at the point {q,p) of the 
phase space. 

The Husimi distribution of a state lip > is defined over the phase space by: 

hii{q,p) = I < gplV-' > P 

and for G L^{1R), we have: 
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\^{x)fdx^ J J\<qp\^>f^ (12) 

To characterize the functions of L^{T*IR) which are {q,p) representations of a state, it is more convenient to 
introduce a complex representation of the phase space z = -^^{q + ip)- Another (proportional) expression of the 
coherent state is then: 

\z >= exp(2;a+)|0 > 

with |0 > being the fundamental of the harmonic oscillator Hq = <f +0^, and a+ being the associated creation 
operator. Indeed: 

,,p>=exp(»--^^)|z> 

The following antiholomorphic function of z is called the Bargmann distribution of ?/;: 

b^{z) =< z\%}) > 

Clearly, we have 

hi,{q,p) |6^(z)|^e~^^, 

hence the zeroes of the function h^{q,p) are those of the holomorphic function 6^(z), which are localized zeroes in 
the phase space. Moreover, ( p^ implies that ij) £ L^{1R) if and only if h^, £ L^{(P, e"'^' /'') and is antiholomorphic. 

The same definitions can be applied for a state >6 Hn{9i,92). The corresponding Bargmann function is a 
theta-function [ p9| and the Husimi distribution is bi-periodic in {q,p), hence is well defined on the Torus Tgp. In this 

representation, Hn is a subspacc of Lf^^{(r,e^^^^ (we keep the weight e"'^''/'' since it is the natural one in the 
Bargmann representation) . 

Remark 2 The Bargmann and Husimi representation of a function of Hn{Oi, 62) are characterized by their N zeroes 
in a given cell Tqp up to a multiplicative constant. Since the N zeroes are constrained to have a fixed sum ( J7^^, we 
get the right dimension — 1 + 1 for the Hilbert space Hn{6i, 62). 



IV. SEMI-CLASSICAL EXPECTATION OF THE CHERN INDEX 



The question is now to guess the value of the Chern index of a given band from classical informations. The first 
result in this direction is a characterization of the Chern index from the Husimi distribution [ p^ : 

Proposition 3 // there exists some point {q,p) £ Tqp of the phase space, such that 

y{di,02) e Te, h^^(^g^^g^){q,p) ^ 

then Cn = 0. 

The proof is quite simple: if < gp|'!/'n(^i, ^2) > is never zero then we can select an eigenstate |V'n(^i,^2) > in each 
fiber such that arg(< qp\ipn{f^i,02) >) = 0, giving us a non-vanishing section of the bundle. This section is also a 
global frame, hence the bundle is trivial: C„ = 0. 

As a corollary, we get an important semi-classical result about Chern indices: 

Proposition 4 // the energy level = {{q,p) such that H{q,p) = E} consists of a single contractible trajectory, 
then the Chern index of the bands of energy around E are semiclassically zero (more precisely, if we consider a h- 
parametrized family of energy bands tending to E as h tends to 0, then for h sufficiently small, the Chern index of 
the energy band must be 0). 

Indeed the WKB construction of quasimodes in phase space ( |30j, [Q) shows that the Husimi distribution 

h^^(^g^ g^-^{q,p) of the eigenstates are highly concentrated and non-vanishing in the vicinity of the trajectory. Thus, 
taking {q,p) on the classical trajectory and from the proposition 2, we obtain C„ = 0. 

Hence, to get a non-zero Chern index, we must investigate situations where the energy level is not connected 
or is non-contractible. In this paper, we will study the first situation. The second one is slightly different, but very 
interesting since it should explain the generic non-zero Chern index arising from Eq. (pT|). 

We need some results about non-zero Chern indices. In fact, for any fixed point {q,p), the Chern index C„ is the 
algebraic number of time that a zero of h^^(^g^^g^){q,p) crosses {q,p) p^ : 
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Proposition 5 Assume that for each 9 G Tg, the eigenvalues En{9i,62) are non- degenerate. Let Z{9i,92) be the set 
of the non-ordered N zeroes of h^^(^g-^ g^y Fix a point {q,p) of the phase space. Define: 

Niq,p) - {(^1,02) e Te I {q,p) e Z(0i,02)} 

Then: 

where the sign ± corresponds to the local orientation of the mapping Zi at {9i,92), where Z = {Zi, Zj^} and 
Z,i9,,92)=^iq,p). 

For example, if the energy level is made of two connected components Fi and F2, then from the tunnelling effect, 
the eigenstates |V'n(6'i, 6*2) > are a superposition of quasimodes > and \1jj2 > localized on each trajectory. If this 
superposition is fluctuating enough when {9i,92) are varying, then h^^(g^fi^){q,p) can vanish for every point {q,p) 
and we expect to get C„ ^ 0. In fact this is possible only for special configurations of the two trajectories. This 
has been investigated in detail in |]lO|, showing that in specific situations, we can observe C„ = ±1. In the following 
sections, we will study the (more complicated) case of three contractible trajectories of energy E. 

V. GENERIC FAMILY OF HAMILTONIAN //(^^ ^^), CHERN INDICES AND DEGENERACIES. 

As in Eq. (|l3|) is not so useful to compute analytically the Chern indices of a given Hamiltonian H. The 
strategy we will adopt is to build a path of Hamiltonians H{t),t £ [0, 1] such that H{0) ~ H and H{\) has zero Chern 
indices (for H{1) no tunnelling effect occurs so eigenfunctions are supported by only one connected and contractible 
trajectory in the phase space, hence Chern indices are zero). Along the path H(t), Chern indices are constant because 
a continuous application with discrete values is constant. Exception is when a degeneracy occurs between eigenvalues. 
In this case the Chern index changes by ±1. In order to calculate the Chern indices, we are therefore left to compute 
these degeneracies. 

In this section, we will study a generic parametrized Hamiltonian family (in other words, we will consider a 
submanifold of the space of Hamiltonians), our main interest is to detect eigenvalues degeneracies. 
There is an essential property for our investigations ( po]]): 

Proposition 6 For a generic family of Hermitian matrices, degeneracies between two levels of the spectrum occur 
with codimension 3. 

This means for example that for a generic 3 dimensional family of Hermitian matrices ^(71,72,73) the value of 
7ii72,73 for which two levels are degenerate, are points in the space of (71,72,73)- 

Consider now a parametrized family of classical Hamiltonian on the torus with external parameters 7 = 
(71, ...,7p). The Fourier components c„i^„2 of H in Eq. (|^), the shape and position of each trajectory in the phase 
space depend on these classical parameters 71, 7p. On the quantum side, the matrix of the Hamiltonian in a specific 
base, depends on the classical parameters 71, 7^ and on the 2 quantum parameters 0i, 02- Since degeneracies are of 
codimension three in the space {9i, ^2,71, 7p), they are of codimension 1 (hypersurfaces) if we project them onto the 
space of classical parameters (71, ...,7^). In fact for each point (71, ...,7p) not belonging to such a hypersurface, the 
Chern index of a given band C„ may be calculated with Eq. (p^. If we cross a degeneracy hypersurface corresponding 
to the band n and n+ 1, the value of C„ and Cn+i changes respectively by ±1 and =Fl (because the sum is conserved: 
cf. [^). For example, in a one- dimensional space 7, degeneracies appear as points. For a two-dimensional space, we 
summarize the previous results as: 

Proposition 7 In a two dimensional space (71,72), degeneracies appear as lines bordering different values of the 
Chern index. The Chern index changes by ±1 when crossing a line. 

In the previous section, we mentioned that these lines occur only if the tunnelling effect occurs. In the next section, 
we will determine (in the semi-classical limit) the location of these lines when the tunnelling effect occurs between 
three trajectories in a periodic cell. 
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VI. DEGENERACIES DUE TO THE TUNNELLING EFFECT BETWEEN THREE TRAJECTORIES OF 

ENERGY E 

In this section, we consider a generic family of periodic Hamiltonians H^. We will use local coordinates 7 = 
(71, ■■•,7p)- 

We assume that for a given value of E, the intersection of the periodic torus with the energy level Y.e consists of 
three contractible trajectories Ti (figure p. 0): 

n T„ = {{q,p) e Tqp/H{q,p) = £;} = Ti u r2 u Tg 

We do not assume that E is non-critical since construction of quasimodes along critical trajectories may be done as 
in H], 0. The only characteristic we will use is that the mean distance between energy levels is of order ln?i| for 
a critical trajectory, whereas it is of order h for a non-critical trajectory. 
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FIG. 2. The energy level E in the phase space. 

Our purpose is now to establish under which conditions (shape and locations of the trajectories Fj) degeneracies 
in the spectrum of H may occur semi-classically near the energy E. More precisely, we will describe the generic 
degeneracy lines in the space 7. 



A. The tunnelling interaction matrix. 

In this section, we construct a basis of Hn{6i, O2) and describe the asymptotic of the matrix of the Hamiltonian H 
in this basis. To have a more precise picture of these asymptotics, we will make assumptions on the respective decays 
of interaction terms. 

For a first reading, one can skip this section and refer to the results mentionned in proposition |^. 

Using the eUipticity oi H — E outside the classical region H{q,p) = E, it is easy to prove that the eigenfunctions 
are of order 0{h^) for all N G IN (noted 0{h°°)) outside the classical region. This property is in fact sufficient to 
construct a basis of i?Ar(0i, 6*2). 

It is easy to modify the Hamiltonian _ff in a new hamiltonian Hi such that: 

• the energy shell E oi Hi is Fi, 

• H = Hi outside the shaded region, see figure H, p.^. 
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FIG. 3. Modification of H. 

Let |1 > be an eigenfunction of li\ corresponding to an eigenvalue E\(K) such that E\(Ji) E a.s h 0. Then 
{H — Ei)\l >= {H — >. In the last expression, the operator is microlocally supported in the shaded region 

where |1 > is not. Hence: 

{H -Ei)\l 0(h°°) 

The same construction applies for the trajectories r2 and and we get the functions |2 > and |3 > microlocally 
supported on r2 and Fa respectively . In the sequel, we will denote by r the greatest of the quantities \ {H — Ei)\i > \ 
(i = 1,2,3). 

^From |1 >, 1 2 > and |3 > we construct Floquet-periodic quasimodes: 

\ipi>= P{ei,92)\i> for i = 1,2,3 (14) 
where P{9i, 6*2) is the operator from L^{R) to Hf^{9i, 62) which makes a given state Floquet-periodic: 

P(0i,02):= eM-^niOi - in202)r^'T^' (15) 

Figure H pM and S pJ| shows the microlocal support of |1 > and \(pi >. 
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FIG. 4. Microlocal support of the states |1 > in the phase space. 
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FIG. 5. Microlocal support of the state \ipi > in the phase space. 

Suppose that we are in the resonance situation: Ei(h) — Ej{h) = o{h,) for i,j = 1,2, 3. Then for e > sufficiently 
small the spectrum of H in the interval / = [Ei{h.) — eh, Ei{h) +eh] is made of 3 eigenvalues Ei{h), E2{K) and E^{K), 
and if we denote by 11/(6*) the corresponding spectral projector of Hpf{9), the spectral space is spanned by: 

\^^>=nj{0)\^,>, 1,2,3. 

Note that for a critical trajectory the width of the interval must be £?i/| ln?i|. 

Let C be the circle of center Ei{h) and radius eh in the complex plane. We represent 11/ and identity as: 



IP- 1 
Ur = — I (H - z)-^ dz, Id^ 



2TTi 



c 



2711 



{El ~ z)-^ dz Id, 



c 



and get: 



J c 



> dz 



1 



= 2— y -{H-z)-\Ei~z)-\H-Ei)\^,> dz 
= 0{r/h^){=0{h°^)) 

since \ \{H — z)~^ \ \ < K/ (eh) and \Ei — z\^^ < K/ {eh) for z G C. In the following text, we replace r//i^ by r, the new 
r function has decay properties similar to the original one (r — 0{h°°) and is probably exponentially decaying). 
Hence > is microlocalized on Ti like \ipi >. 

The next step is to get an orthonormalized basis from the Qipi >)i. Let 



Since N is positive hermitian, there exists a unique positive hermitian matrix B such that: 



Let B = (with f3ij = f3ji) and: 



B-^ = N 



\ei >=^l3,j\'ipj > 



(16) 
(17) 

(18) 



We claim that the family je; > is orthonormal. Indeed: 

< ei\ej > = ^l3ik(3ji < > 

k,l 

= {BNB),, 
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The tunnelling interacting matrix is the 3x3 Hermitian matrix: 

A={<e,\H\e,>),, (19) 

The eigenvalues of A are exactly the eigenvalues of H which belong to the / interval. 
We will now compute the semi-classical asymptotics of Oy . Let: 

G= (<^,|i/l^j >)„=i,2,3 (20) 

^From Eq. (p^), we get: 

A = BGB (21) 

hence we want to compute asymptotics of B and G. 

First we remark that in Eq.(|l6|), we can replace \'ip > by \(p >, since by Pythagoras theorem: 

< (p^\ipj >=< ■4},\i}.j > + < v?j - Vil'/'j - >=< > +0{r^) (22) 

Since H commutes with 11/, we get by the same method: 

< ipt\H\fj >=< ipi\H\tljj > + <if,- ipz\H\ipj - tpj >=< ipilHltpj > +0{r'^) 

We have 

< ipi\ipj >= 6ij + 0{r) 

Hence: 

N = ld + 0{r), B = ld + 0{r) (23) 

Now, we want to expand the functions \ipi > using ( p^ ) and ( [Tsl ) and evaluate the scalar products. ^From Eq.(|20|) 
and (|lj), we get: 

< Vil-fflV-i >=< MH\1 > +2^Re (exp(-m^) < l\H\ln >) (24) 

where: 

• n = (ni,n2), = (6*1, 6*2), 

• D = {{Z xIN)\ {~N X {0})} is the half plan of n, 

• \lfi >— Tq^Tp^ll > is a quasi mode concentrated on the translated trajectory {Ti)ni,n2 in the (^1,712) cell of 
the phase space. 

The first term of ( |2^ ) is Ei + 0{r^). The term < l|i/|lf{ > is of order 0{r) because it corresponds to the 
tunnelling interaction between the quasimode |1 > localized in cell (0,0) and the quasi mode \lfi > localized in the 
cell (?T,i,n2) 7^ (0,0) (we may need to modify the function r to get this decay: r denotes the strongest tunneling 
interaction between two components of the energy shell). Hence: 

< ^i\H\iJi >= Ei+0{r) 

Similar asymptotics hold for the other diagonal terms. 
A non-diagonal term of G is e.g.: 

< Vi|-ff|V'2 >= exp(-m^) < l\H\2,-i > (25) 

And there is generically only one leading term due to the strongest tunnelling interaction between |1 > and 12^5^2 > 
(located in the cell ni2: ni2 may be for example (0,0) or (0, ±1) or (±1,0)), of strength 

< l|i?|2Hj2 >= mi2exp(icpi2) 
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where mi2 > is 0(r). This gives: 

< >= mi2 exp(i(v3i2 - ^126')) 

and similar expressions for others non-diagonal terms. Hence all non-diagonal terms are 0{r). If we denote by f the 
second strongest tunneling interaction, we have 

< il)i\H\ip2 mi2exp(i((pi2 - 71126')) + 0(f) 

We assume of course that f — 0{h°°)r. 

^From A — BGB, B = N^^/'^ and N = Id + 0{r), we get for the diagonal terms of A: 

au ^ E, + 0{r) 

and for the non-diagonal terms: 

ttij = ruij exp(i(</3y - fiijO)) + 0{f) (26) 
To get a more accurate asymptotic of the diagonal terms, we will now study the shifted interaction matrix: 

< e,\H ~ E^le, B{G - EzN)B 

We have G — E3N — (< ipi\H — E^lipj >)ij. For the non-diagonal terms, we get similar asymptotics as above ( |2^ ) 
with new constants. For simplicity, we will however keep the same notations, since the decays remain of the same 
order and the phases were not explicit. For the diagonal terms, we get by the same method: 

(< tp.i\H ~ Esltp, >)^E^-E3+ mu cos{(pu - nuO) + ... (27) 

Since B — I modulo 0{r) and G — E^N is of order 0(r), the shifted interaction matrix is G — E3N modulo an error 
of order 0{f). 

In this section we have used he fact that the tunnelling interaction gives terms of order r — 0{h°^). In fact the 
reader may replace 0{h°°) by 0(e^'^/'*) which is conjectured to be the right decay order, where S is similar to the 
Agmon distance between two trajectories (as in the Schrodinger case H — + V{x)). See for partial 

results in this direction). 

We summarize the results of this section: 

Proposition 8 Let H be a Hamiltonian and E an energy such that the energy shell H = E is made of three contractible 
connected components. Let |1>, |2>, |3> be the quasimodes of H localized near one of the three classical trajectories 
corresponding to eigenvalues Ei{h) — * E, E2{K) —> E, E:i{h) — > E Suppose that Ei{h) — Ej{h) = o{K). Then for £ > 
and h> sufficently small, the spectrum of H acting on Hn{9i, 62) in the interval [i?3(?i) — eh, E3{h) + eh] is made 
of three eigenvalues. 

We have constructed an orthonormal basis je^ > spanning the corresponding 3-dimensionnal vector space. Let r be 
a majorant of the largest tunneling interaction between two different trajectories, and r be a majorant of the second 
largest tunneling interaction between two wells (r = 0{h°°) and f = 0{h°°r). More precisely, we choose i and j 
m[l,3], we look at the second strongest interaction between trajectory i and translated of the trajectories j, excluding 
the translation if i — j . r and r are conjectured to be exponentially decreasing with respect to r). Then the matrix 
of H in the basis \ei > has the following asymptotic: 

• The non- diagonal terms are 0{r), more precisely: 

= ruij expi^ipij - rtij.B) + 0(f), 0(r) 

• The diagonal terms are given by: 

an = Ei + ma cos{ipii - n^.O) + 0{f) 
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If we look at the dependency of the interaction matrix with respect to the external parameters 7, it is easy to prove 
that: 

d^Oij = {d^Ei)5.ij + 0{r) 

where dj denote the differential. The first term is dominant, since the variation of Ei is of order h (or h/\ ln?i| for a 
critical trajectory). 

Since we are interested in eigenvalue degeneracies, we may substract -E1./3 to the interaction matrix, hence the 
variation of the shifted interaction matrix with respect to 7 is described by the variation of Ei — E3 and E2 — E3 with 
respect to 7 up to an error of order r. Hence, we will reduce our parameter space to two parameters 71 and 72, and 
we justify in appendix ^ that they can be chosen as: 

71 = (^1 - E3)/h (28) 

72 - {E2 - E3)/n 

(replace fi by T%/\ ln?i| for a critical trajectory) . 

Moreover, since the error is of order 0{r), the degeneracy lines in the parameter space will be of 0{h °") s ize (and 
conjectured to be of exponentially small size). We will come back to this point more precisely in section \VI C . 



B. Model of a 3x3 interaction matrix, and computation of its Chern indices. 

This paragraph is self-contained. From the previous paragraph, we have to consider a continuous mapping on the 
torus Tg = [0, 27rp into the 3x3 Hermitian matrices: 



where the diagonal terms are: 
and the non-diagonal terms are 



A : (01,92) eTg^ A{0i,92) = (a.,.j)»,,=i,2,3 

a = Ei + ma cos{ipii - fiiiO) i^ 1,2,3 (29) 

j — mij expi(i^ij — fiijO) i < j — 1,2,3 (30) 



with fiij = (?T.i(ij),n2(ij)) £ and — {0i,02)- We consider (m^ , (^.^ , ny )ij=i.2,3 as fixed. We will study only the 
dependance of A with the parameters {Ei,E2, E3). 

We recall that (niij , tpij , fiij) refer to the tunnelling interaction between the three trajectories {Ti)i,i = 1,2,3 as 
sketched in figure (|^), and that Ei is the energy of the quasimodes on trajectory F^. 

If no degeneracy occurs in the spectrum of ^(6*1, 6*2) (for all = {0i,02) G Tg), each eigenvector family has a well 
defined Chern index. Precisely, each eigenvector family {\ipn{0) >)g for n = 1, 2, 3 is a submanifold of the projective 
space P{(r^), homeomorphic to Tg. The complex line bundle structure of P((Z'^) induces a complex line bundle over this 
submanifold whose topology is characterized by its Chern index. For specific values of the parameters (codimension 
1 in the space of {Ei)i) degeneracies occur and this causes a change of the Chern index. In this paragraph we will 
compute these Chern indices and the locus of the degeneracies in the space of {Ei)i for each family of eigenvector 
{\4'n{0) >)g for n = 1,2,3. First we remark that substracting the diagonal matrix E^Id to A does not change its 
eigenvectors, so the only non-trivial external parameters are: 

71 = (^1 - E3)/n (31) 

72 = {E2 - E3)/h 

The locus of the degeneracies we are looking for are lines in the space (71,72)- 

A general property of 3 x 3 Hermitian matrices, mentioned in the appendix ^ is that a degeneracy occurs for the 
matrix A if and only if: 

9(012023031) = (32) 
ail - 0-22 H = 

ai3 0,23 

, 013032 012023 

a22 - 033 H = 

ai2 ai3 
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9 stands for imaginary part. Moreover, the degeneracy is between levels n = 1,2 (respect, n = 2, 3) if K(ai2a23a32) > 
(respect < 0). 

The first equation can be written: 



e arg(ai2a23a3i) = ipu + (p23 + Vsi - (^12 + ^23 + n3i)9 = [n] 



(33) 



The phase 8 can be seen as the total tunnelling phase for the cycle of trajectories (ri,r2,r3). We want to find 
solutions in 9 of this equation for fixed values of (pij, fiij. For that purpose, we assume that {ni2 + ^23 + n^i) 7^ 0. 
This is a generic assumption. This means that the cycle of trajectories (ri,r2,r3) is not contractible on the torus 
T 

To simplify notations, we will assume that ni2 = ^23 — 0, ^31 — (0,1) and na — (1,0). This corresponds to the 
tunnelling interactions sketched with dashed lines in figure (^. For the tunnelling problem considered in this paper, 
we think (without proof) that this case is the general one after a suitable lattice transformation in SL{2, Z). But for 
the self-contained problem of this paragraph, the general case is solved in the appendix |^. 

P 



® 



FIG. 6. Example of tunnelling interaction. ni2 — n2-i ~ 0, nai = (0, 1) and nn — n22 ~ W33 ~ (1, 0). 



Eq. ( p3| ) has then two solutions 02 = <y3i2 + ^23 + fai [tt] 
Let 



s = exp(i9) = ±1. 



The two-last equations of Eq. (02) give: 



El - E3 = smis ( — !^ - ) + 77133 COs((/)33 - 6*1) - TTln cos((^ii - 61) 

\m23 mi2/ 



(34) 



and similarly for E2 ^ E^- 
Let us remark that: 



A cos((?i — a) — B cos(0i — /3) = (A cos a — B cos j3) cos di + {A sin a — B sin (3) sin 9i 

= D cos(6'i — 77) 



where: 



So we put: 



(cosry, sin 77) = 



D = yjA^ + 52 - 2AB cos(a - /?), 

A cos a — B cos /3 A sin a — B sin /? 



(35) 



D 



D 



Ki — \/ mf-^ + 77133 ~ 2miim33 cos tpi 



(36) 



K2 — V 7^22 + 1T^33 ~ 2m22m33 COS iy92 



where: 
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and define: 



il = TOl3 i2 = TO23 

\m23 mi2/ VTO13 mi2/ 



From and (^^, i/ie degeneracy lines in the space (71,72) are </ie following t E M parametrized curves (s = H): 

jiit)^{Kicos{t^7j)-sYi)/h (37) 
j2it)^{K2Cos{t)-sY2)/h 

where s — ztl, t = di — ri2. rj = rj^ — ri2 and the angles 7/2 cmd 773 may be computed using (p^. 

As t describes [0,27r], (71 (t), 72(0) describes two translated ellipses of axes parallel to 71 = ±72, one for each 
s e {— L 1}- s = 1 (respect, s = —1) gives the degeneracy line between levels n = 1, 2 (respect, levels n = 2, 3). See 
figure (0). These two ellipses may intersect, but they don't really intersect in the whole space (6*1,6*2,71,72) because 
they correspond to two different values of 62- 

Outside the ellipses, the Chern indices are zero. This is because, when 71,72 — > 00 the A matrix goes to a diagonal 
matrix with trivial eigenvectors. Crossing a degeneracy line changes a Chern index change by ±1. It can be shown 
that the value ±1 is related to the orientation relatively to the orientation of the parametrized ellipse. From this, we 
deduce the values of the Chern indices in figure (|^) . 






n=l n=2 n=3 

FIG. 7. Degeneracies lines and value of Chern indices from Eq.(^7|). 



C. The degeneracy lines are ellipses. 

We are now ready to present the final result. Consider a generic two-dimensional submanifold of the manifold 
of classical Hamiltonian H-y (7 are coordinates on this manifold, and can be viewed as external parameters of the 
Hamiltonian) . We have obtained that in this space, degeneracy lines are two ellipses (one for each pair of levels (1, 2) 
and (2,3)). 

Note that it make sense to state that the shape of the curves are ellipse even if the i?^-space is a manifold. This is 
because the ellipse are exponentially small, and thus live in the tangent space. 

More precisely, there exist a coordinate system (71,72) (i.e. a special parametrization of the Hamiltonian) for which 
the two ellipses are given by Eq. (^7|). The width of the ellipses are given by the coefhcients if 1,^2, whereas the 
position of the centers are given by coefficients sYi, SY2, s = ±1. 

It is obvious from Eq. (pfl) that Ki, K2 are exponentially small in the semi-classical limit. But this is true also for 
Fi , 12 • This is shown at the beginning of the appendix (^ . 

We have then two different generic situations in the semi-classical limit: 
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case: the non-diagonal terms 1^,12 are exponentially small with respect to the diagonal terms Ki,K2- The 
two ellipses are identical, see figure || p. At the origin, the three bands of energy will have Chern indices of 
(=F1, ±2, =f1)- Outside the ellipse, the Chern indices are (0,0,0). Note that in the (71,72) parameter space, the 
ellipse are exponentially small and exponentially flat because Ki is exponentially small compare to K2 (or the 
inverse) . 




n=l 






v'2 





/ > 







n=2 n=3 
FIG. 8. Degeneracies lines and value of Chern indices in case 1. 



the diagonal terms Ki, K2 are exponentially small with respect to the non-diagonal ones Yi, 12- See figure 
At the scale Yi of the distance between the two center, the degeneracy lines collapse to 2 single points 
Individually these lines are still elliptical with respect to their own scale Ki. 



case 

Bp._ 

as fi tends to 
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6^ 



n=l 







^0 



(9' 



n=2 n=3 
FIG. 9. Degeneracies lines and value of Chern indices in case 2. 



The non-generic intermediate case occurs when the diagonal terms Ki and non-diagonal terms Yi are of the same 
order. Hence the degeneracies are two ellipses (s = ±1) with symmetric centers. These two ellipses may intersect. 
Corrections to the leading behaviour in the semi-classical limit can modify the shape and sizes of the ellipses, see the 
numerical section. 

We now discuss the special translation-symmetric case corresponding to the first calculation of Chern indices done 
by Thouless and al. [0. For example, take H{q,p) from Eq. ( ^9|) with 6 = 71 = 72 = 0. See also remark]^. We 
construct |2 >= Tp/all > and |3 >= rp/3|2 > by translation by P/3 in the momentum direction. Applying the 



relation Tp/^Tq 



e^^-^I^TQTpi^ we get: 



(/?2 = </'i/2 



2-kN 



(38) 



Since the non-diagonal m^j are equal, so we are in the first case (double ellipse centered at the origin of figure | 
Assuming that 3 doesn't divide N then we get K\ = K2 and the precise shape of the ellipse: 



K 



71 



cos(t), 72 



-cos(i±-) 



We conclude that we get an ellipse of axes 71 = ±72 and excentricity \/2/3. 

Remark 9 We want to stress that all the results obtained in this section are valid modulo error terms of order 0{r^), 
where r = 0{h°°) is the greatest tunnelling interaction between two wells. If one of the previous contributions mij 
between two other wells is smaller than r"^ , it is not significant. For example, if we consider an Hamiltonian where 
wells 2 and 3 are very closed compared to well 1, our method must be modified, although we think that our results Eq. 
are still qualitatively correct. In this situation, a more convenient approach is to first treat the cluster of wells 2 
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and 3 and finally consider the interaction between the cluster and well 1. By this approach, we obtain results similar 
to those described in JT^/ for the 2-wells problem. 



VII. NUMERICAL ILLUSTRATION 



More precise informations on how to perform the numerical si mulation, and a software will be ava ilable on the web 



address: http : //www-f ourier . uj f -grenoble . f r/ ~parisse or ittp : //lpm2c . polycnrs-gre . f r/ 



f aure 



We take the following Hamiltonian on the torus, parametrized by 71 and 72: 

H{q,p) = cos(27r^) + cos(67rf ) 

+ 45sin(27r^)cos(27r|r) 
+ ((71 - 72) cos(2^^) + ((-71 - 72) sin(2^^) 



(39) 



The two terms in the first line are the main part of the Hamiltonian. They create three symmetric wells in the p 
direction in the lower part of the energy spectrum, see figure (10). The second line breaks this symmetry, by moving 
the second well on the left when b increases, b can be seen as a perturbation from the translational symmetric case. 
The coefficients 71 and 72 in the third line change the depth of wells 1 and 3 respectively. We therefore expect them 
to act directely on (£'1 — E2) and {E^ — E2) respectively, such that (Bl) will be satisfied. In figure ( p^ it seems clear 
that the main tunnelling interactions between the wells are those sketched in figure (^) as assumed for the calculations 
of the previous section. 



FIG. 10. Contour levels of Hamiltonian (691), with h = 0.2 and 71 



-72 



-0.1 



For the numerical calculations, the Hamiltonian (39) has been diagonalized and formula ( |l0|) has been used to 
obtain the Chern indices C„. Maps of the Chern indices have been numerically obtained by varying parameters 71 
and 72. 

For iV = 17 the values of the Chern indices for the first three energy bands n = 1, 2, 3 are shown in figure |Tl|. 
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n=l 



n=2 



n=3 




2 2 




FIG. 11. Chern indices for the first tfiree bands n — 1, 2, 3, for A'^ — 17. Case (a): The strength of the perturbation is fo = 0.1. 
Case (b): The strength of the perturbation is 6 = 0.2. 

In case (a), the strength of the perturbation is b — 0.1. The two elhpses are degenerated. This is the generic 
situation depicted in figure and the same generic situation as for the translation-symmetric case with b = 0, 
treated in the previous section. This means that the parameter b = 0.1 is small enough to stay in the same generic 
ensemble. 

In case (b), the strength of the perturbation is & = 0.2. The two ellipses are separated, as in the second generic 
situation, sketched in figure (^. 

In case (a) and (b), the two ellipses have the same shape, in agreement with the semi-classical results of the previous 
section. 

Figure |l^ shows results for a lower value of A^, = 11, and for strengths b — 0.2 and b = 0.25. Here, the generic 
situations expected in the semi-classical limit N — > +oo are less visible. There are still elliptical curves (from a 
topological point of view), but their position and shapes are affected by non-leading corrections with respect to 
N +00. 
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-o!l5 -0.1 -Q.I -q!i5 -0.1 



FIG. 12. Chern indices for the first tfiree bands n = 1,2, 3, for A'^ = 11. Case (a): The strength of the perturbation is 6 = 0.2. 
Case (b): The strength of the perturbation is b = 0.25. 

VIII. CONCLUSION 

In this paper, we have shown in which conditions a non-zero Quantum Hall conductivity can occur in the framework 
of the Harper model, and for the tunnelling between three trajectories in a periodic cell. In this framework, the Hall 
conductivity is proportional to the topological Chern index. These results were derived semi-classically and give at 
the same time a description of degeneracies in the spectrum. 

Currently, one tries to observe experimental signatures of the Harper spectrum (Landau level substructures) in 
lateral superlattices with periods of about lOOnm on GaAs-AlGaAs hetero junctions [[l2|j24| ]. Our results could therefore 
have some experimental importance in the future. 

The main mathematical gap in this work is that we do not have results about microlocal tunnelling effect (like 
those obtained by Helffer and Sjostrand ||l3| in the Schrodinger case). The main ideas to solve these difficulties are 
currently niajorations of the wave function (A.Martinez [pl|), complex paths (estimations of the tunnelling effect by 
Wilkinson |]3l| , ^ ), or normal forms (S.Nakamura [^), but the authors are not aware of a general rigorous result. 

Some other related problems could be investigated in the future like: 

• numerical localization in the 7-space of the ellipsis (this requires numeric estimations of the tunnelling effect) , 

• generic interaction between N bands (A^-wells tunnelling effect) by a recursive clustering approach, 

• the tunnelling effect between non-contractible classical trajectories 
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APPENDIX A: EIGENVALUE DEGENERACIES IN 3x3 HERMITIAN MATRICES. 

In the text, we use the following result of linear algebra that we have not found in the literature: 
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Theorem 10 Let ^ &e a 3 x 3 Hermitian matrix: 



A 



a S J 
She 
7 e c 



where a,b,c are real numbers and 7,(5, e are non-zero complex numbers. 
Then A has a eigenvalue A with multiplicity at least 2 if and only if: 

a ~ b =c = A e R 

e -f 

(this is a system of Z real equations) . The third eigenvalue is X — tr{A) — 2A and the spectrum is (A, A, A') with A' > A 
z#3?(fc7) >0. 

Proof 

Let us recall that a hermitian matrix has real eigenvalues and is diagonalizable. Then A has an eigenvalue A with 
multiplicity 2 or 3 if and only if its minimal polynomial has degree 2 or 1. The matrix is then a linear combination 
of / and A. We have: 

/a^ + \S\^ + h\^ (a + 6)5 + 76 {a + c)-f + Se 
A^ = i * 62 + |£|2 + |(5|2 {b + c)e + h 

V * * |7|' + |e|' + c2 

where the stars denote the conjugate of the symmetric coefficients. 

The other eigenvalue is A' = tr{A) — 2A, where tr{A) denotes the trace of the matrix A. Therefore the minimal 
polynomial is: 

P{X) = {X -X){X - tr(A) + 2A) =X^ + {X- tr(A))X - (2A - tr(A))A 
A^ = (tr(A) - X)A + A(2A - tr{A)) 



hence: 



Eventually we get the system: 



f a2 + + |7|2 = {triA)-X)a + X{2X - tr A) 



{tr{A)~X)b + X{2X-trA) 



a' + |£p + |7|2 = (tr(A)-A)c + X{2X - tr A) 

{a + b)5 + 7£ = (tr(A) - X)5 

(a + c)7 + 5e = {tr{A) - A)7 

L {b + c)e + 6"/ = (tr(A) - A)e 



(Al) 



But tr{A) — a + b + c, so the last three equations simplify to: 

■'ye — (c — X)6 
6e = (6-A)7 
Sj = (a — A)e 

When S, 7, and e are non-zero, we verify that these three last equations imply the first three of system Eq. ( |Al| ). 
Indeed by multiplying two with two, we obtain: 

|ep = (&-A)(c-A) 
\5\' = ia-X){b-X) . 
|7p = (a-A)(c-A) 



The spectrum is (A, A, A') iff 3A < tr{A). This gives 3?((5e7) > 0. 
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APPENDIX B: INFLUENCE OF THE EXTERNAL PARAMETERS 



In this appendix^ we show that genericaUy, only two external parameters 71 , 72 control degeneracies and the Chern 
indices, as in Eq. (28). 

From the generic cyclic assumption made in section VI B| , we have the property: 

7^ * 7^ k,j 7^ fc, TnijTTijk is 0{h°°) with respect to rriik 

In the analogous Schrodinger situation, this means that the triangular inequality is strict for the Agmon distance. 

Remark 11 For the Agmon distance, the converse assumption is also generic (wells may be "shaded" by other wells). 
It is excluded here, because assuming that a well is shaded implies that the cycle of trajectories is contractible. 

For example, if H is symmetric under translation Tp/^, we see that: 

m.ijmjk/mik = m^- 0{r) = 0{h°°) 

Then we get easily that Ei - Ej must be 0{h°°) in Eq. (||). 

More precisely, we want to look at the 7-dependency of the two equations of (^4|). From the first equation, we 
define 



f "^12 _ 



smi3 + 77133 cos(v333 - 6*1) - mil cos(i^ii - 6*1) 

\TO23 mi2 



/l(7,^l) = ^((^l-^3) 

and we define /2 similarly using the second equation of so that (|3j) becomes 

/(7,^i) = 

where 



/ 

We see easily that: 



RP X [0, 2tt] 

(1,0) ^ (/l,/2) 



df = UiEi - E3,E2 - E3) + 0{h^) 

h 

If we want to have the simplest possible parametrization of the degeneracies, we have to find 2 parameters such that: 

= ^2 parameters (^(-^1 ~ ^3), -{E2 - E3)) 

is invertible. Hence, we can not take 61 as a parameter, and we must take two 7 parameters. This means that the 
degeneracy will not be described by a point on a 1 parameter-7 line but by a line on a 2 parameters-7 plane. We 
now assume that the parameters 71 and 72 satisfy the invertibilty hypothesis, and more precisely that there exists a 
constant C such that: 

|J|<C, \J-^\<C, withJ = a^„^,(l(^l-^3),^(^2-^3)) (BI) 

for 7 in an neighboorhood of a point 7*^. In addition, if we suppose that 7*^ is an approximate solution of / = 0: 

/(7^ei) = 0(/i°°). 

we can apply the implicit function theorem with parameter ?i, 0i, 73, 7p, for h £]0, Hq] with Hq sufficiently small such 
that the rest terms 0{h^) arc smaller than the constants C). We get the existence of an ?i-independent neighboorhoud 
U of ((7*')i, (7*')2) such that for every fixed (73, 7p) in a neighboorhoud of ((7°)3, (7°)p) and for every fixed 9i, 
the function /(.,73, ...,7j,,6'i) of (71,72) is bijective from U to f{U). Since /(7°,6'i) — 0{h°°), there exists a unique 
point 7(6*1,73, ...,7p) in U such that /(7(6'i,73, •■■,7p),^i) = and !7(6'i,73, ...,7p) - 7"] = 0{h°^). 

If we fix (73, ....,7p) and let 9i move, the projection of the point 7(^1,73, •.•,7p) in the (71,72) space describes a 
curve of size 0{h°°). At that scale, all classical quantities (like tunneling interactions) are constant (with a relative 
error of order 0(/i°°)). 
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APPENDIX C: MORE GENERAL DEGENERACY CURVES 



In this appendix, we solve the general problem raised in section VI B, i.e. we look for the expression of the degeneracy 
curves in the space X2, X3 for any value of n^j, i,j — 1, 2, 3. 
Let us put 



The first equation of Eq. (|32|) gives: 



e = (f- Ne = Q 



(CI) 



with ip = (pi2 + (p23 + ^31- We put s = exp(i6). We suppose (see paragraph ( |VI B[ ) ) that 

N ^6 

There is an unique decomposition of iV as iV = pAI where p £ N* is the greatest common divisor of (iVi, A^2), so 
(Mi,M2) are relatively prime. Then Eq. (CI) gives 



M9 = i(p + kn) /p, 



,2p 



(C2) 



k even (respect, odd) corresponds to s = +1 and degeneracy between n = 1,2 (respect, s = —1 and degeneracy 
between n = 2, 3). 

The last two equations of Eq. (B2) give the curves equations as: 



71 = h{-sYi + mil cos(v3ii - nii9) - 771,33 cos((^33 - 71^36')) 

72 = h{-sY2 + 17122 C0S(V322 - "-226*) ^ "I33 COs((^33 - 7133^')) 



(C3) 
(C4) 



with 71, Yi defined in paragraph ( VIB ) and 9 is allowed to vary but constrained by Eq. (|C2|). Because M = (Mi, M2) 
are relatively prime, there exists no from Bezout 's theorem such that M and ng form a basis of the Z'^ lattice, that 
is det[M, 7?o] = 1. In this basis, we can decompose: ftu — biM + CifiQ, with bi, Ci G Z for i = 1, 2, 3. This gives 



with 



hi^p/p - bik2TT/p, k = l,..,2p 



1,2,3 



and a free parameter t — no9 G iR.The final expression for the degeneracy curves is: 



71 (i) = fi{-sYi + mil cos(ri^fc - cit) - 777,33 cos(T3,fc - czt)) 

72 (i) = fl{-sY2 + 77722 COs(r2^A; - C2t) ~ 777,33 C0s(r3,fc - Czt)) 



(C5) 
(C6) 



with /s = 1, ..,2p. This gives 2p curves. They can be quite complicated in general. See an example in figure (|13|). 
These curves are degenerate if all ci — Q, this means that every 77^^ is proportionnal to N . As pointed out in paragraph 
( |VI B[) , we think that the only possible situation in the Harper model considered in this paper is the ellipse curve, 
Eq.(|37D, for which p = 1, ci = C2 = C3. 
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